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We present a general model study of surface-enhanced resonant Raman scattering and fluores- 
cence, focusing on the interplay between electromagnetic (EM) effects and the molecular dynamics as 
treated by a density matrix calculation. The model molecule has two electronic levels, is affected by 
radiative and non-radiative damping mechanisms, and a Franck-Condon mechanism yields electron- 
vibration coupling. The coupling between the molecule and the electromagnetic fleld is enhanced 
by placing it between two Ag nanoparticles. The results show that the Raman scattering cross 
section can, for realistic parameter values, increase by some 10 orders of magnitude (to ~ lO"'^* 
cm'^) compared with the free-space case. Also the fluorescence cross section grows with increasing 
EM enhancement, however, at a slower rate, and this increase eventually stalls when non-radiative 
decay processes become important. Finally, we flnd that anti-Stokes Raman scattering is possible 
with strong incident laser intensities, ^ 1 mW//im. 

PACS numbers: 33.20.Fb, 78.67.-n 33.50.-j, 42.50.-p 



I. INTRODUCTION 

Surface-enhanced Raman scattering (SERS) was dis- 
covered three decades ago. SERS attracted a lot of at- 
tention through the mid-1980's>ii2i^ During this period 
the basic mechanisms behind the effect were studied, 
explained, and debated. The advent of single-molecule 
(SM) SERS*'^'?^'! started a second period of intense re- 
search. Now probably the prospect of utilizing SERS and 
related spectroscopic techniques as an extremely sensi- 
tive analytic tool, possibly combined with scanning probe 
techniques)^ in a variety of life-science applications^ pro- 
vides the main motivation for research in the field. But 
SM SERS experiments performed with intense lasers has 
also raised new questions about the fundamental mecha- 
nisms involvediiSiiiiiSiH 

It is generally agreed that electromagnetic (EM) en- 
hancement effects are the most important reason for the 
dramatic increase of the Raman scattering cross-section 
GR seen in SERS experimentsi^J^iiS*iLiS*i2i2a In addi- 
tion to this, may also be enhanced due to charge 
transfer effects iSLSSiSiSi The EM enhancement effects 
received a lot of attention from the theory side in the 
early days of SERS and has continued to do so until 
now. The electromagnetic enhancement has in the gen- 
eral case not one single reason, but involves several, more 
or less closely related aspects such as plasmon resonances, 
lightning rod effects, and the formation of "hot-spots" in 
fractal clusters. 

However, single-molecule Raman scattering is not pos- 
sible without a molecule that scatters the laser light 
inelastically. In this paper we focus on the interplay 
between the EM enhancement and molecular dynam- 



ics, a topic that has received relatively limited atten- 
tion in the literature so far. A brief account of this 
work was published in Ref. We treat the molecule 
as an electronic two-level system. Thus, we make no at- 
tempt of calculating ab initio molecular properties, and 
charge transfer processes are also outside the scope of 
this work. The focus is instead on calculating how the 
molecule's different states are populated and how its co- 
herent dipole moment develops given a certain energy- 
level structure, electron-vibration coupling, electromag- 
netic enhancement and laser intensity. By treating the 
molecule dynamics within a density-matrix calculation, 
we can evaluate a combined fluorescence and Raman 
spectrum where also effects of electromagnetic enhance- 
ment are directly included. Density-matrix methods were 
employed by Shen?® to distinguish between Raman scat- 
tering and hot luminescence, but to the best of our knowl- 
edge they have not been used in the context of SERS. 

The results show how, for a molecule placed between 
two metallic nanoparticles, both the fluorescence and 
in particular the Raman cross sections are much larger 
than for a molecule in free space. This can be dis- 
cussed in terms of two EM enhancement factors, M and 
\Md\. Given an electric-field enhancement M(a;) at the 
position of the molecule, the Raman cross section in- 
creases by a factor ~ \M\'^ compared with in free space, 
whereas the fluorescence cross section increases by the 
factor ~ |Af |''/|Mdp. |Mdp is a measure of how much 
the decay rate of an excited state of the molecule is am- 
plified near one or several metal particles. For moderate 
to large molecule-particle distances M and \Md\ are al- 
most equal, but when the molecule gets very close to 
a particle (a few nm or less) \Md\ can be much larger 
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than M . Thus, surface enhancement of fluorescence is 
much less marked than that of Raman scattering, and for 
small enough molecule-particle separations the fluores- 
cence cross section saturates. When studying a molecule 
next to a single metal particle, M does not reach at all 
as high values as in the two-particle case, but \Md\ still 
does. Consequently, the SERS effect is largely absent in 
this case, and fluorescence is strongly quenched when the 
molecule is placed close to the particle. Distancing the 
molecule from the particle, \Mj\ drops and the fluores- 
cence goes through a maximum and then eventually falls 
back to the free-molecule value. 

We have also studied how the Raman spectrum de- 
velops when the intensity of the driving laser field is in- 
creased to relatively large values of order lmW/(/im^). 
In this case we find that anti-Stokes Raman scattering 
becomes possible even if the molecule vibrations are not 
thermally excited. Such effects have been observed in 
a number of SERS experiments, but the issue has been 
quite controversialii2iii*i^ We find that anti-Stokes Ra- 
man scattering in our model occurs because a vibra- 
tionally excited level in the electronic ground state is pop- 
ulated by (possibly repeated) excitation from the laser 
and subsequent deexcitation. We also find that when the 
laser intensity is increased to values where Rabi oscil- 
lations become important, the anti-Stokes Raman cross 
section decreases, and the corresponding peak in the 
spectrum broadens. 

The rest of the paper is organized as follows. In Sec.HTI 
we describe the model for the molecule that we use. As 
already stated, the electromagnetic enhancement plays 
an important role in SERS, and Sec. IIIII describes how 
the EM enhancement of the incident laser field, the emit- 
ted light, and the deexcitation of the molecule is calcu- 
lated. Moreover, we also present numerical results for 
some representative cases. Then in Sec. |^ we go back 
to the model molecule and investigate how the parame- 
ter values entering the model affect the absorption and 
Raman scattering cross sections for the molecule in free 
space. In Sec. we use density-matrix theory to derive a 
general expression from which the combined fluorescence 
and Raman scattering spectrum can be calculated for a 
molecule placed near metal nanoparticles, and driven by, 
in principle, an arbitrarily strong laser field. The so cal- 
culated spectra are presented in Sec. IVII and Sec. IVIII 
concludes the paper with a brief comparison with ex- 
periments. Finally, Appendix IXI gives some background 
information about the electromagnetic calculations, and 
Appendix ^ presents a derivation of the non-radiative 
damping of the molecule due to electron-hole pair cre- 
ation in the metal particles. 



II. MODEL 

A schematic illustration of our model, both in terms of 
molecule-nanoparticle geometry and the main features of 
the molecule model is given in Fig. ^ In this figure we 
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FIG. 1: Schematic illustration of the model we use. (a) Elec- 
tromagnetic enhancement is achieved by (usually) placing the 
model molecule between two metallic nanoparticles. (b) The 
molecule has two electronic states (sometimes called bands in 
the following) each with a number of vibrational sublevels. 
The equilibrium position for the vibrational coordinate is dis- 
placed by xo upon electronic excitation; this mechanism pro- 
vides electron- vibration coupling. The arrows indicate typical 
pathways for the molecule state in fluorescence and Raman 
processes. 



also indicate the main steps involved in different Raman 
and fluorescence processes. 

We model the electronic degrees of freedom of the 
molecule as a two- level system with states \g) and |e). 
The frequency fi^e , deflned from the energy- level separa- 
tion — Eg as 

nrige^E.-Eg, (1) 

and the dipole matrix element, expressed as the prod- 
uct of the elementary charge Cg and a dipole length £dip, 
Po = egidip, (to be further discussed below) are the only 
parameters we need in this context to characterize the 
essential electronic properties of the molecule. 

To deal with Raman scattering we must of course in- 
troduce vibrational degrees of freedom and an electron- 
vibration coupling. We consider only one symmetric, vi- 
brational mode, and let Q denote the corresponding co- 
ordinate. The vibrational mode is characterized by its 
angular frequency fivib and reduced mass /i. The purely 
vibrational Hamiltonian can thus be written 

ffvib = 7i^^vib(6^fo+^), (2) 

where b and are annihilation and creation operators 
for a vibrational quantum. We assume that the vibra- 
tional frequency is independent of the electron state, but 
the equilibrium position is displaced a distance Xq in the 
excited state. Moreover, the transition dipole moment 
introduced above is actually a function of Q: in the Born- 
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Oppenheimer approximation we get the Taylor expansion and 



p EE p{Q) = {e{Q)\e,z\g{Q)) ^ po + Q + . . . , 

around Q = 0. (Thus, we focus on the z components of 
the electric fields and dipole moments, since they are the 
ones enhanced in the nanoparticle geometry we study.) 
Then the dipole matrix element between two electron- 
vibration product states yields 

{xo]m\{e\eez\g)0;n) 



{Q-n\{g\H'\e)\x^-m) = -(poSo/2)e'"^* /(71, m). (9) 

(In the following, when the local field at the molecule is 
modified by the nearby metallic particles, these matrix 
elements must be corrected by an enhancement factor.) 

The interaction between the molecule and the vacuum 
fluctuations of the electromagnetic field can be described 
by a Hamiltonian -fffiuct = —SezEvac ■ z, where the corre- 
sponding electric field must be given on second-quantized 
form. In free space we have 



:po(a:o;m|0;n)+ (^^) (xq; m|Q|0; n) -f . . . , (3) t) ^ ^ ^ZAI^^^^ 

k.a 



where \xo',n) denotes vibrational state n in an oscilla- 
tor potential with the equilibrium position displaced to 
Q = xo, etc. We adopt the Condon approximation, thus 
retaining only the first (Albrecht A) term,^^ in Eq. 
It gives the dominating contribution in situations where 
the incident light is close to resonance with the electronic 
transition. We thus get 

{xq; m\ {e\eez\g)Q; n) = paf{n, m), (4) 

where the Franck-Condon factor 

f{n,m)= {0;n\xQ;m) = {xo;m\0;n) 



(_l)(™ + fe)Q," + ™- 

^ k\{n ^ ky.{m ^ k)\ 



(5) 



and the ratio a between xo and the average zero-point 
vibration, (((2Q)2))i/2 ^ {2n/{fin^iy,)y/^ 



a = Xq 



2h '■ 



(6) 



serves as a measure of the electron-vibration coupling in 
the model. 

The molecule is furthermore interacting with the 
electromagnetic field, both with the incident laser field 
and the electromagnetic vacuum fluctuations that cause 
spontaneous emission. The electric field originating from 
the laser can be written 



El = zEoCosVlLt = z{EqI2) [e 



(7) 



so that the nominal incident intensity is /in = ceQE'^/2. 
We model the laser- molecule interaction within the dipole 
approximation by a — Ce r- E term, H' = — eei?o-z[e*^'"* -I- 
e^'^^*]/2 in the Hamiltonian. When evaluating matrix 
elements due to H' we adopt the rotating wave approx- 
imation (RWA). Only the part of the field varying with 
time as e~*^^* is kept in the matrix elements when the 
molecule is excited, and vice versa when it is deexcited. 
This yields 

(xo;m|(e|i/'|g)|0;n) = -(poi?o/2)e-'"-*/(n,m), (8) 



(10) 



where the electromagnetic field has been normalized in a 
box with side L, and ak,a and ^ are annihilation and 
creation operators for photons. For a molecule in free 
space the interaction with the vacuum fluctuations gives 
a transition rate 



7l■ad,o(?^,™) 



37r?ieoC' 



■bon/(n,m)|^ 



(11) 



where lo is the angular frequency corresponding to the 
transition energy, i.e. 



u ~ Qge + (m — n)wvib- 



(12) 



With a dipole moment corresponding to £dip = 1 A, a 
transition energy of 2.5 eV, and |/p ~ 1, we get numer- 
ically 7rad,0 ~ 5.9 X 10^ s"^ 

This transition rate is enhanced when the molecule is 
placed near metallic nanoparticles; the radiative losses 
increase and in addition energy can be dissipated in the 
particles, thus 

7rad,o("-, »7l) |Af<j(tj)p7rad,o(?T-,TO)- 

We will discuss the dissipation enhancement factor 
|Mrf(tj)|2 further in Sec. Hill 

Our model also includes, at a phenomenological level, 
a few more relaxation mechanisms. Vibrational damp- 
ing is described by the parameters 7^f,^ and 7^^^. The 
transition rate from a state with n vibrational quanta to 



the 



ith 



1 is given by wy^fl and n'yl'rl 



b' 



the 



electronic ground and excited states, respectively. The 
electronic state of the molecule does not change in this 
process. Finally, we include a dephasing rate jph affect- 
ing the coherent dipole moment of the molecule. The 
primary effect of the dephasing rate on the calculated 
results is to broaden the fluorescence resonances of the 
molecule. In reality an organic molecule has many vi- 
bration modes, and therefore an almost continuous fluo- 
rescence spectrum. Dephasing gives us a broadened flu- 
orescence spectrum even though the model molecule has 
only one vibrational mode. Similar broadening param- 
eters are also used in Raman cross section calculations 
that account for molecular structure in more detailj^ 
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III. ELECTROMAGNETIC ENHANCEMENT 
A. Theory 

In this section we address the calculation of the electro- 
magnetic enhancement factors M and \Md\ introduced 
above. We first assume that the system of metal par- 
ticles and the molecule is illuminated by a p polarized 
plane wave with angular frequency co and wave number 
k = uj/c incident from a direction specified by the angles 
9 and (p. The corresponding incident electric field can be 
written 



and 



(13) 



with 



Eq — Eo{x coa9 coaip + y cos 9 sin ip — zsin6'), (14) 

(in the following we assume the time-dependence of all 
fields to be e"*"^*). When this field impinges on the metal 
particles they are polarized, plasmons may be excited, 
and, most importantly, the electromagnetic interaction 
between the spheres will cause a field enhancement in the 
region of space in between them. We are mainly inter- 
ested in calculating the z component of the electric field 
at the position of the molecule (r = 0) ; it is this field that 
excites the dipole moment of the molecule when a laser 
beam illuminates the system. We define the enhance- 
ment factor M (uj) in terms of the induced total field at 
the position of the molecule, through the relation 



M{uj,9)^z-Etot{6)/\Eo\. 



(15) 



The enhancement of the incident laser field is thus given 
by AI{ilL,6in)- Moreover, as a consequence of electro- 
magnetic reciprocity, the field sent out in the direction 
of 9 and ip by an oscillating dipole (angular frequency lo) 
placed at r = is equally enhanced by a factor M{u!,9). 

To calculate the local field and M we must find the 
electric field around the nanoparticles. To this end 
we employ extended Mie theory, expanding the field 
around each of the two spheres in terms of vector 
spherical harmonics representing magnetic and electric 
multipoleS )^^i^° 



Eir) 



E 



to) 



rs) 



(16) 

Here the functions V'*-^'' and ?/'^'^-' represent waves that 
are regular at the origin and outgoing from the sphere, 
respectively. The index s tells whether we are referring 
to the upper (s = 1) sphere, centered at fi, or the lower 
one (s = 2). The rest of the indices indicate the values 
of the angular momentum of the multipole £ and m, and 
whether it is a magnetic (r = 1) or electric multipole 
(r = 2). The corresponding basis functions are given in 
terms of vector spherical harmonics Ximf 
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'ipiim = Xim{9,p)zi{kr) 



(17) 



V^2£r„ = fc-iy X |x,™(0,(^)zi(fcr)} , (18) 



where zi is either a spherical Bessel function ji for the 
regular waves, or a spherical Hankel function hg, for the 
outgoing waves. The expression in Eq. (|15|l is valid for 
points situated between sphere s and a larger concentric 
sphere that just touches the other sphere. For a general 
point outside the two spheres, on the other hand, the field 
must be written as a sum of the outgoing waves generated 
by both the spheres, plus the incident wave driving the 
system. 

The a and b coefficients appearing in Eq. (|16|l are re- 
lated by sphere response functions su and S2£ depending 
on the radius R of the sphere and its dielectric proper- 
ties characterized by the local dielectric function e (taken 
from Ref. ISlf) and wave number kj. — y^'euj/c. We get 



bu 
an 



kRj[{kR)^ji{kR)Ji 



kRh\{kR) - hi{kR)Ji' 

for the magnetic multipoles, and 

^ 62^ e kRj[{kR) + ji{kR) je-l-Ji) 
a2i> ekRh[{kR) + hi{kR){e-l-Ji)' 

for the electric multipoles. J7i is shorthand for 



(19) 



(20) 



Jl = krRj[{krR)/il{krR). 

In view of these relations between incident and outgo- 
ing waves, we have full knowledge of the electromagnetic 
field once we know the a coefficients on both spheres. We 
solve for them by realizing that the waves incident on a 
sphere either originate from the external source or from 
the waves scattered off the other sphere, and this yields 



rim 



r^m, cxt 



m.Tim 



(21) 

Here s stands for the other sphere (i.e. 1 = 2 etc.). Ex- 
pressions for the c and aoxt coefficients can be found in 
Appendix ^ 

After solving the system of equations obtained from 
(|21|l . also the b coefficients can be determined thanks to 
Eqs. lfTO|l and (|2nj). Then, as discussed above, the electric 
field can be calculated anywhere in space, in particular 
the local field at the molecule is given by 



(22) 



from which we can find the enhancement factor M . It 
should be noted that in the calculation of the z compo- 
nent of the electric field on the symmetry axis there are 
only contributions from m — Q terms. 

When calculating |Md|, we cannot use the plane wave 
as a source, instead we place an oscillating point dipole 
at the position of the molecule. 



p(t) = Re \pe 



(23) 
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In free space, it would send out radiation with a total 
(time-averaged) power^ 



4 ? 

47r£n Sc-^ ' 



(24) 



With the spheres present we can once again solve for the 
total electric field from Eq. H21|l . The only difference is 
that now Ccxt appropriate for a localized dipole source 
has to be used, see Appendix IXI Once the field has been 
calculated we integrate the time-averaged Poynting vec- 
tor over a small sphere enclosing the dipole to find the 
total radiated power Pmio, which yields a contribution to 
the damping enhancement rate l-M^p, 



|Mrf(co)|2 =PMic/Pf 



free 



(25) 



Pmic accounts for losses of energy due to radiation leaving 
the molecule in all directions, as well as for dielectric 
losses in the particles, i.e. energy from the molecule that 
goes to heating the metallic particles. 

In practice we need to add an another contribution 
to |Mdp that is due to electron- hole pair creation in 
the metal particles when the molecule is placed very 
close to themjSSi^^iM This involves processes that depend 
on the non-local optical response of the metal particles 
and which are therefore not included in the calculation 
discussed above. The dielectric losses captured by the 
Mie calculation at short molecule-particle separations d 
scale as 1/d? as a result of the distance dependence of 
the molecule's dipole field. The power loss caused by 
electron-hole pair creation, on the other hand, scales as 
and becomes a dominant damping process for small 
d {d less than 1 nm or so). In Appendix IbI we give a de- 
tailed presentation of the calculation of the power loss 
Peh due to electron-hole pair creation in the metal par- 
ticles. After Peh is calculated the total damping rate 
enhancement is found as 



\MdiL0)\^ = (PMie + Peh)/Pbee. 



B. Results for the enhancement 



(26) 



Figures 13 El and 0] show calculated results related to 
the electromagnetic enhancement. In these calculations 
the molecular response which we address in more de- 
tail later plays no role, however, the placement of the 
molecule, i.e. where to evaluate the electric fields, is cru- 
cial. We consider a case where the molecule is placed 
symmetrically between the two metal nanoparticles, d is 
the molecule-particle separation, and consequently the 
smallest distance between the spheres is 2c?. 

Figure El shows results for M{uj, 9 = 90°) as a function 
of photon energy for a few different values of d. There is 
an overall increase of the enhancement when the spheres 
approach each other. Moreover, the enhancement factor 
has a resonance peak with a resonance frequency that 
shifts with d. For d = 25 A it lies somewhat above 2.6 eV, 
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FIG. 2: (color online) Calculated enhancement factors M{uj) 
for the incident light as a function of photon energy, for a 
molecule symmetrically placed between two Ag spheres with 
400 A radius. The different curves give results for a number 
of different molecule-particle separations d. 
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FIG. 3: (color online) Calculated enhancement factors |M|^ 
(full curves) and as a function of the molecule-particle 

separation d. As in Fig. |21 the molecule is symmetrically 
placed between two Ag spheres of radius 400 A, and the pho- 
ton energy is 2.2 eV. The dashed curves give |Mdp without 
the contribution due to electron-hole pair creation, whereas 



the dash-dotted curves show 
included. 



Aldl with these contributions 



but it redshifts with decreasing d and for d — 5 A ends up 
at « 2.2 eV. The resonance is caused by the interaction 
between the plasmon modes of the two separate spheres 
leading to the formation of a coupled mode with surface 
charges of opposite sign facing each other across the gap 
between the spheres. This means that a smaller d gives a 
more charge-neutral mode, and thereby smaller restoring 
forces leading to the redshift shown in Fig. [21 

In Fig. O we have plotted the enhancement factors 
|M(a;)p and |M<j(w)p (for the latter quantity both with 
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FIG. 4: Calculated damping rates as a function of photon en- 
ergy. As in Figs.|21and|3tlie molecule is symmetrically placed 
between two 400 A radius silver spheres. The curves show re- 
sults for three different molecule-particle separations d for the 
separate contributions from, on one hand, electron-hole pair 
creation and, on the other hand, losses captured by the Mie 
calculations (here labeled as radiation and ohmic damping). 
The damping rate for the free molecule was calculated from 
Eq. inil with ^dip = 1 A. 



electron-hole pair losses constitute a minor contribution 
for d = 20 A, but has at d = 5A become the dominant 
damping mechanism. 



IV. THE MOLECULE IN FREE SPACE 

Given the model for the molecule discussed above, we 
can calculate the absorption and Raman scattering cross- 
sections for a molecule in free space using lowest-order 
perturbation theory, i.e. the Fermi golden rule»2& 

The absorption cross-section derived in this way reads. 



2 \ ^ 



?irtot/(0,n)/(0,n) 



''■ ^ tot 



(27) 

where Ftot = Iph + 7vib/2 + 7rad,o/2 and 7rad,o and 7vib 
are the radiative and vibrational decay rates of the fi- 
nal, excited state. For the combinations of parameter 
values that we use, jph gives the completely dominating 
contribution to Ftot for a molecule in free space. 

For the (fundamental) Raman cross-section we get in 
a similar way, using the Fermi golden rule with a second- 
order transition matrix element, 



and without contributions due to electron-hole pair cre- 
ation) as a function of the molecule-particle separation 
d for photon energy huj = 2.2 eV. Quite naturally, since 
d spans a large range of values, so do the enhancement 
factors, from about 1 for d > 1000 A, to 10^ or more 
at the smallest d < 10 A. For distances d larger than 
200 A the three curves follow each other very closely; at 
least on this scale no difference is discernible. This is be- 
cause the silver particles are close enough to enhance the 
incident field and thereby, as a consequence of electro- 
magnetic reciprocity, also enhance the radiation rate and 
\Md\- But the particles are still sufficiently far away from 
each other and the molecule that the dielectric losses are 
negligible. Thus, radiation losses completely dominate 
other damping mechanisms here. Continuing towards 
smaller d, the two iM^p curves separate themselves from 
the |Mp curve because now losses in the silver particles 
are no longer negligible compared with radiation losses. 
The inclusion of damping due to electron-hole pair cre- 
ation does not make any difference at first; non-local ef- 
fects play no role and all loss mechanisms are already 
accounted for in the Mie calculation. It is only when d 
reaches values of « 30 A or below that the two curves 
representing lAf^jp with and without electron- hole pair 
damping begin to differ. 

Figure0]shows the separate contributions to the damp- 
ing rate from electron-hole pair creation and the remain- 
ing, radiative and dielectric loss mechanisms. As in the 
previous figure we see a rapid increase of the damping 
rate with decreasing d. This tendency is more marked 
for the electron- hole pair losses which vary as The 



c^i? = Po 



OJ 



2^4 



E 



/(l,^)/(0,n) 



ihuj 



vib 



(28) 

Here uj' = Ql - Wvib, and Ftot = Jph + 7vib/2 + 7rad,o/2 
for the intermediate, virtual state. 

Turning to calculated results for spectra, we begin by 
looking at absorption and Raman profiles and their de- 
pendence on the molecular parameter values. Figure [51 
shows absorption and Raman profiles plotted for a series 
of different dephasing rates jph- In panel (a) we also give 
experimental results for the absorption cross section of 
a commonly used fluorescent dye molecule, Rhodamine 
6G (R6G). This comparison provides some indications on 
what parameter values can be considered realistic. Thus 
we have chosen Tiflge — 2.35 eV to get about the same 
peak position in the model calculation as for R6G. The 
value used for hflvih, 160 meV, is characteristic for a 
C-C stretch vibration. The dipole length £dip and the 
dephasing rate largely determine the height and width 
of the absorption peak. We have set £dip = 1.2 A, while 
we in the calculations presented in Fig.jSjused a number 
of different values for •jp/^ to study its effects. Obviously 
7ph — 1-3 X 10"'^^ s~^ gives the best agreement with the 
R6G experimental spectrum; the remaining difference is 
the broader wings (especially on the low-frequency side) 
of the calculated spectrum. Finally, the Franck-Condon 
parameter a mainly determines the strength of absorp- 
tion sidebands (or shoulders) due to vibrational excita- 
tions. In the calculations presented here we have used 
a = 0.5 and this gives a shoulder in the absorption spec- 
trum of similar strength to that found in the R6G spec- 
trum. 
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FIG. 5: (color online) Calculated (a) absorption and (b) Ra- 
man profiles (i.e. cross sections as a function of incident pho- 
ton energy) for a molecule in free space for a number of dif- 
ferent dephasing rates 7ph. The remaining parameter values 



used in the calculation are £dip 
Tifivib = 160 meV, and a = 0.5. 



1.2 A, ft^ge = 2.35 eV, 



As can be seen in Fig.|Sl(a), varying results in vari- 
ations of the width and height of the absorption peaks; 
the widths are proportional to jph whereas the heights 
are approximately inversely proportional to jph- Thus, 
with 7pft considerably smaller than 1.3 x 10^^ s~^ the 
vibrational sidebands create marked, resonant peaks at 
2.51 eV and 2.67 eV, while with a larger jph the absorp- 
tion spectrum is more broadened. 

The Raman profiles displayed in Fig. O (b) show the 
same trends as the absorption profiles when the dephas- 
ing rate jph is varied, but the dependence is somewhat 
more complicated in this case. The Raman scattering 
cross section is governed by a second-order matrix ele- 
ment with an energy denominator with a real part set 



by the laser detuning and an imaginary part mainly set 
by jph- Thus in the range of laser frequencies close to 
resonance we get peaks in the profile with a height deter- 
mined by 7p/i. In this situation [in Fig. [S] (b) primarily 
between 2.3 eV and 2.5 eV] the height varies as ~ '^/iph 
and the width as ~ 7^^. On the other hand, further 
away from resonance the detuning, not the dephasing 
rate, dominates the energy denominator entering the Ra- 
man scattering matrix element and in this case an varies 
much more weakly with 7p^. 



V. CROSS-SECTION CALCULATION 

In this section we outline in detail the calculation of the 
spectrum of light emitted by the molecule. This includes 
both light scattered as a result of Rayleigh or Raman 
processes, and fluorescence as a result of electronic tran- 
sitions in the molecule. The methods presented below 
makes it possible to carry out this calculations in a gen- 
eral case with strong enhancement of both the incident 
light and of the damping rate. 



A. Emitted intensity 

The emitted light intensity at the point ro in the far 
field can be writtcn'^^ 

Jen, = 2eoc{E'f\vo.t)E^+\vo,t)). (29) 

Here -E'g^'' (ro ^ t) and Eg (tq , t) stand for the positive and 
negative frequency parts of the 9 component of the elec- 
tric field, respectively. By using the normal-ordered cor- 
relation function {Eg \ro,t) Eg^\rQ,t)) we are assured 
that the vacuum fluctuations do not contribute to /cm- 

Our present goal is not only to find the total inten- 
sity of emitted light, but also its spectral distribution. 
Starting from Eq. (|29|l . and using the Wiener-Khintchine 
theorem this can be written 

1 

= -Re / dT2eoc{El-\ro,0)El+\ro,T))e''^\ 
Jo 

(30) 

The electric fields are caused by the electric dipole mo- 
ment p(t) of the molecule. A point dipole with dipole 
moment zpe~*"* placed at the origin in free space gener- 



ates the electric field E = Ope^^'^'^ 



ot) 



'sin 0/(47reoC^r) 



at a distance r from the dipole. If this expression is com- 
bined with Eq. (|30|l we get the differential scattering and 
fiuorescence cross-section 



w"* sin^ I 



Re 



dre'""(p(-)(0)p(+)(r)). 

(31) 



When the molecule is no longer placed in free space 
the electromagnetic propagation from source to detector 
is modified. This can be described by the enhancement 
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factor M{uj,9) introduced in Sec. IIIII M{uj,9) is closely 
related to a photon Green's function.'^^ The sin 6 factor in 
Eq. H31|l is absorbed in M which yields the emitted spec- 
trum (in the numerical results presented below we use 
6 = 90° as the observation angle) from a molecule placed 
in the vicinity of one or several metal nanoparticles, 



d^a _ uj*\M{uj,e)\^ 



■Re / dre'"^^ {p^-~\0)/+Ht)). 

(32) 

To get any further we must evaluate the dipole-dipole 
correlation function {p^~^ (0) p^"*"' (t)) by solving the equa- 
tions of motion for the molecule dynamics. We will deal 
with this using density-matrix methods. 



different vibrational levels are included in Ctr in a phe- 
nomenological way. The last term CphP makes it possible 
to introduce additional phase relaxation, also in a phe- 
nomenological way. 

The molecular Hamiltonian is diagonal, and can be 
written as a sum of electronic and vibrational energies, 

A^vib-l 

Hmoi^ J2 ^ \l;n> (Ei+nhn^ih) <l;n\. (34) 

To see more clearly what H' in Eq. H33|l means, let us 
focus on a manageable example with just 4 levels (iVvib = 
2). We then have 



B. Equation of motion for the density matrix 

We consider a model of the molecule with in total N 
quantum states: two electronic levels, \g) and |e), some- 
times called bands in the following, with iVvib vibrational 
states per electronic level, i.e. N = 2Nvih- The density 
matrix p is thus a 2A'vib x 2A'vib matrix with equation of 
motion 



H',p] + CtrP + CphP- 



(33) 



The first term to the right governs the motion as a re- 
sult of the molecule Hamiltonian i?moi, and the interac- 
tion H', between the molecule and the laser field. The 
term CtrP yields the damping of the density matrix as 
a result of transitions caused by iJfluct in which photons 
are spontaneously emitted, but also transitions between 



v+/(o,o) y+/(i,o) 
14/(0,1) y+/(i,i) 



l/_/(0,0) y-/(o,i) 
V.f{l,0) y-/(i,i) 



and H' is explicitly time-dependent since 
14 = -poEoe^'''^' /2. 



(35) 



(36) 



With these equations, and their generalizations to cases 
with more vibrational levels, the first term in Eq. I|33() 
can be calculated. 

Consider now the second term £trP, in Eq. ()33|1 . Given 
a spontaneous transition rate r^j (at zero temperature) 
from level j to level fc, (in our model this rate is due to 
vibrational damping for intraband transitions and due 
to radiative damping and electron-hole pair creation for 
interband transitions), this term can be written^ 



J 



CtrP = - y^(»fcj 
jk 



PO-jkO-kj - 2cFkjP(Tjk\ - y^(n 



r 



]k 



kj 



[(^kjO'jkP + pcrk](^jk - ^CTjkPcrkj] 



(37) 



Here aab stands for an operator (or matrix) with the ab 
element equal to 1, and all the other elements equal to 
0. These operators fulfill the relations 

(^ab(^kl = (^alhki and [(yab,(^kl\ — (^alhk — CrkbSal- (38) 

The first sum in Eq. (|37|l refers to (possibly thermally 
activated) decay, the second to thermal excitation, and 



(rifej) = l/[e 



{E,-Ek)/{kBT) 



(39) 



implying that the thermal bath consists of a number of 
harmonic oscillators, something that is at least certainly 
true for the radiation damping. Given that all excitation 
energies Ej — Ek occurring in our model are fairly large 
compared with fc^T, we have restricted the calculations 
to the zero-temperature limit, (rikj) = 0. 



The last term in the equation of motion describes de- 
phasing of interband (i.e. ground-excited and excited- 
ground) coherences. This only involves one density- 
matrix element at a time, thus 



Cph P — Cph ^ O-ij Pij — — «7p/i ^(^ij Pij + <^ji Pji)j 
ij ieg 

(40) 

where i G g indicates that when the molecule is in state 
i it should be in the electronic ground state, etc. 
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C. Stationary state 



D. Calculation of the spectrum 



We need to solve Eq. and begin by determining 

the stationary-state populations and coherences. To fa- 
cilitate the solution in a case with an arbitrary number 
of levels we form an A'^^-dimensional vector p, from the 
elements of the density matrix as 

P= [pii,P2i, ■ • ■ Pni,Pi2, ■ ■ ■ Pnn]- (41) 
The equation of motion can then be written 

• I^Za (42) 

where the tensor £ describes the coupling between two 
density matrix elements caused by the Hamiltonian and 
the damping. C can be deduced from the right hand side 
of Eq. (ESI- 

The stationary-state density matrix pss is essentially 
time-independent, but because of the explicitly time- 
dependent terms in the Hamiltonian, it does have co- 
herences (dipole moments) oscillating with the laser fre- 
quency. We therefore make the ansatz 



Pss 



e Po, 



(43) 



where pq and are time-independent, and is a diagonal 
tensor. Moreover, the diagonal elements of fl referring to 
populations or intraband coherences equal 0. Only the 
diagonal elements of referring to interband coherences 
are non-zero. They equal -t-^L, for "excited-ground" co- 
herences and —flL for "ground-excited" coherences. For 
a two-state system (just electronic degrees of freedom) 
we would have 



n = 



- 








0" 
























. 








0. 



(44) 



The stationary state can thus be determined by solving 
the equation 



. dpss 
'' dt 



Cpss ^ e-'^"r!po = /:e-^^"po, (45) 



which after a multiplication by e'^* can be rewritten as 
[£' - n]pa = 0, (46) 

where 

£' = e*^*£e-*^* (47) 

is time-independent. (Basically, we have transformed the 
equation of motion to an interaction representation.) In 
addition to Eq. (|46|l . po niust satisfy Tr[po] = 1- 



We can now focus on calculating the key quantities: 
the fluorescence and Rayleigh and Raman scattering 
cross-sections, as given by Eq. The positive fre- 

quency part of the dipole operator originates from tran- 
sitions from an electronically excited state to the elec- 
tronic ground state. The negative frequency part, on the 
other hand, comes from transitions from the electronic 
ground state to the excited state. With the aid of the a 
operators introduced before Eq. H38() we can then write 
the expectation value appearing in Eq. (|32|l . 



aj'Ge 
|2 ^ ^ 

b,keg 



(p(-)(0)p(+)(r)) = bop E fib,a)fik,j){aabiOWjiT)). 

(48) 

The expectation value {aab{0)akj{T)) can be evaluated 
using the quantum regression theorem (QRT), as we will 
show nexti^ 

With the QRT the value of any element of the density 
matrix at time r can be expressed as a linear combina- 
tion of the density matrix elements at the earlier time 0, 
provided that the process is Markovian. For a physical 
process to be Markovian there must not be any memory 
effects, damping should be frequency- independent, and 
the perturbing noise completely "white." It is fairly clear 
that none of these conditions are fulfilled in a strict sense 
here, for resonant enhancement of the coupling between 
the molecule and the electromagnetic field implies that 
damping and noise is stronger at some frequencies than at 
others. However, the EM enhancement, see Fig. 12 varies 
relatively slowly, on an energy scale AE of a few tenths 
of cV, corresponding to a memory time scale h/AE of 
a few femtoseconds, shorter than the other time scales 
relevant for the molecule dynamics. At the same time we 
note that in a number of interesting problems in quan- 
tum optics, for example an atom in a photonic crystal 
with an atomic transition energy near a band-gap edge 
of the photonic crystal, involves very important memory 
effects that require other theoretical methods and leads 
to novel physical phenomenai^Si^ 

Continuing the calculation, we first consider a one- 
operator expectation value and note that 



W> = Trb(T)crfcj] = Pjkir). 



(49) 



We express this dependence through a Green's function 

Gjk,rs{T), thus 



P]k[T) = E Gjfe,rs(T)prs(0). 



(50) 



In the previously used vector and tensor notation this 
equation would read p{t) = G{t)p{0). In terms of oper- 
ator expectation values Eq. H50I) means 



(^fc.M) =E^^■fc.-(^)('^-(0))■ 



(51) 
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The quantum regression theorem2& states that the re- 
lation between the two-operator correlation functions 
(o'a6(0)crfcj(T)) and {(Tab{0)<7sr{0)) IS identical to the one 
between one-operator expectation values expressed by 
Eq. l5ll), thus 

{<Jab{0)<Jkj{T)) =^Gjfc,,,(r)((7a6(0)(7,,(0)). (52) 

rs 

It will now be possible to calculate the dipole-moment 
correlation function since the expectation value appear- 
ing in Eq. (|52(l can be evaluated from the stationary-state 
density matrix pss, and the Green's function, in view of 
Eq. (|5()|l , can be deduced from the equation of motion for 
the density matrix. 

To explicitly calculate the Green's functions we need to 
solve the equations of motion for the density matrix given 
certain initial conditions at time t = 0, since obviously 
Eq. states that Gji,ki{T) is the value that the ele- 
ment pji will take at time r given that all density matrix 
elements are at time t = except the pki which equals 
1. [This should be viewed strictly mathematically; the 
initial conditions here sometimes mean that the density 
matrix is traceless and non-Hermitian. However, when G 
operates on a physical p{0) one obtains a sensible result 
also for p(t).] We can thus write down an equation of 
motion for the Green's tensor 



dG 
(— - 

dt 



CG^ilS{t), 



(53) 



where the right hand side takes care of the initial condi- 
tions [i.e. G(t) = for i < 0]. 

To solve Eq. (|53|l we make an ansatz analogous to the 
one made for pss above. 



Git) 



-iCtt 



Go{t). 



(54) 



Then by introducing the Fourier transform Go{uj) of 
Go(t) through 



Go(c^) = / dte'-*Go(t), 

we arrive at the solution 

Ga{ijj) ^ + i5)i + Vl -2'] \ 



(55) 



(56) 



From a formal point of view the imaginary part 5 — Q, 
but to avoid that the solution diverges at the driving 
frequency $7^ we introduce a finite, but small, 5. This 
gives a width to the calculated Rayleigh scattering peak. 
In the calculations reported in this paper we have used 
the value 5 = 2 x 10^^ s""'^ unless otherwise stated. 

We can now insert the expression for the expectation 
value into the equation for the cross-section, 

oj^\M{io,e)\\ ,2"^%. ,,,, 

-\Po\ 2^ fib,a)f{k,j) 



b,keg 



/>oo 

;Re / dTe'^-y^G,k.,rs{r){(Jab{0)'Jsrm. (57) 



The time integral gives the Green's function Fourier 
transform. In view of Eq. H54|l and the fact that the 
index j refers to an excited state while k refers to the 
electronic ground state the time integral can be written 



dr e*""e-*^'^"Gojfc,.«(r) = GojkM^ " ^l). (58) 
Furthermore, the expectation value can be simplified to 

(crafc(O)crsr(O)) = {(yar{0))Sbs = PQ,raSbs, (59) 

which finally yields 

uj^\M{lo,9)\\ ,2"^%. ,,,, 
- bo! 2^ f{b,a)f{k,j) 



b,k£g 



X PO,ra^'^[Gojk,rb{^ " ^l)] (60) 

r 

for the scattering and fluorescence cross-section. 

E. Results in an elementary case 

Equation Ht)0|) contains contributions to the cross sec- 
tion due to both diagonal and off-diagonal elements of the 
density matrix. These contributions typically have very 
different physical origins. The light emission related to 
diagonal elements are due to an excited state being pop- 
ulated before decaying; this describes fluorescence pro- 
cesses. The off-diagonal elements, on the other hand, 
represent an oscillating dipole moment on the molecule, 
and the corresponding contributions to the emitted light 
are due to various scattering processes. 

To illustrate this we show results for the scattering 
cross section in a very simple case in Fig. |H| We con- 
sider a model molecule with only two electronic levels 
without any sublevel structure due to vibrations (cf. Ref. 
40). The electronic excitation energy is set to 2.3 eV and 
the incident laser photon energy is 2.33 eV. The elec- 
tromagnetic enhancement is used as a parameter: both 
|Af|^ and \Md\'^ are set to 10^, independent of the fre- 
quency u. As for damping mechanisms, we of course 
keep the radiation damping, but there is no dephasing in 
this model. By inspecting Eq. H60I) we see that indices b, 
k, a, and j in this case refer to one definite state (whereas 
they in the general case run over a set of different vibra- 
tional states), only the index r can point to either the 
ground state or the excited state, so the cross section is 
a sum of two terms. The three curves in Fig. El show 
the total cross section as well as the contributions from 
the term proportional to the element po,ee (diagonal) and 
the term proportional to po.ge (off-diagonal). The total 
cross-section is everywhere positive as it should, and has 
a sharp peak at htu = 2.33 eV due to Rayleigh scattering, 
(the width of this peak is set by the parameter S used in 
the calculation) . The Rayleigh scattering is of course the 
result of photons being re-emitted because the molecule 
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FIG. 6: (color online) The spectrum of scattered light from 
a two-level system. In addition to showing the total cross 
section with a Raleigh peak as calculated from Eq. 161)1 . the 
other two curves show the contributions originating from the 
diagonal and off-diagonal elements of the density matrix, re- 
spectively. 



FIG. 7: (color online) Calculated fluorescence and Raman 
spectra for a molecule placed symmetrically d = d\ = d-i be- 
tween two spherical (R = 400 A) Ag nanoparticles. Different 
curves correspond to different particle-molecule separations 
d. The molecule parameter values used are MIl ~ 2.45 eV, 
^itJvib = 160 meV, jph = 1.3 x 10^* s-\ -f^^^^^,) = 2 x lO'" s-\ 
7vib(«) 



a = 0.5, and /i„ « 0.13 nW /{firn^ 



has an oscillating dipole moment described by the off- 
diagonal elements of the density matrix, and this curve 
follows the one for the full cross-section closely near the 
laser frequency. However, for photon energies near the 
electron transition energy the off-diagonal contribution 
becomes negative, largely canceling the diagonal contri- 
butions in this case. 



VI. CALCULATED SPECTRA 

We begin by looking at spectra calculated for a 
molecule placed symmetrically between two spherical Ag 
nanoparticles. Figure [7| displays a number of combined 
fluorescence and Raman spectra corresponding to differ- 
ent molecule-particle separations d, and thereby different 
electromagnetic enhancements. The plotted differential 
cross section corresponds to a situation where both the 
incident and scattered light propagate in the symmetry 
plane (i.e. 9 = 90°). All the spectra show a broad fluores- 
cence background. For the smaller values of d, sharp Ra- 
man peaks rise above this background. The fundamental 
Stokes peaks, red-shifted by TiWvib from the laser photon 
energy Tl^lL, are the highest ones, but at least for d = 5 
A we can also see an overtone peak at hflL — ZTiWvib due 
to creation of multiple vibrational excitations. The very 
strong dependence of the Raman peaks on the molecule- 
particle separation d is the most striking tendency seen 
in the plot. Increasing d from 5 A to 8 A reduces the 
peak height by approximately a factor of 6, at d = 12 A 
only a small Raman peak remains, and at d = 16 A it 
is nearly impossible to discern a Raman peak. The main 



reason for this is the stronger EM enhancement one gets 
with a smaller d, cf. Figs. [3 and|31 As was shown in Fig. 
3 of Ref. l25| the Raman scattering cross section behaves 
as CTfl, « |M(rii)p|i\f (il — Wvib) per i?,, free as long as d is 
not too small. The Raman cross section scales with the 
fourth power of the enhancement because Raman scat- 
tering involves two steps: in the flrst step a photon is 
temporarily absorbed and the molecule goes into a vir- 
tual state, in the second step a photon is emitted while 
the molecule goes back to the ground state, albeit to a 
different vibrational state. The rate of both these steps 
are enhanced by a factor |Af p. 

Also the fluorescence background in Fig. [7| changes 
with d; the fluorescence cross section shows an increasing 
tendency with decreasing d, however, this change is not 
at all as marked as for the Raman signal. Naively one 
may think that fluorescence, which involves an absorp- 
tion event and an emission event, should also display a 
cross section scaling as \M\'^. But this is not so, because 
in the case of fluorescence the molecule is in a real, elec- 
tronically excited state after having absorbed a photon 
but before emitting the fluorescence photon. Looking 
at the final step of the fluorescence process it is then 
clear that there are two factors that determine the fluo- 
rescence intensity: (i) the EM enhancement for the fre- 
quency of the emitted photon |Af (a;)p, and (ii) the prob- 
ability of finding the molecule in the excited electronic 
state Pcxc- The probability Pexc is affected by the electro- 
magnetic enhancement, but in two competing ways that 
largely tend to cancel each other. Poxc increases when the 
molecule is excited from the electronic ground state and 
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FIG. 8: (color online) Calculated spectra for three different 
frequencies of the driving field. The molecule-particle separa- 
tion is d = 5 A, and the rest of the parameter values are the 
same as in Fig.|7| 



the rate of such processes scales as |Af ($1^)^. However, 
at the same time the excited state is emptied by radia- 
tive (including fluorescence) and non-radiative processes, 
and the rate of these is given by |Mdp7o,rad^exc- Con- 
sequently, the probability of finding the molecule in an 
electronically excited state depends in a stationary-state 
situation on the EM enhancement factors as 

Pexc - \M\y\Md\^ 

and for the fluorescence cross section we get 
ap ^ iMpFcxc ~ |Af|V|Afd|2. 

In Fig. |S1 we compare three spectra calculated with 
different laser photon energies. These spectra result 
from the combined effect of a frequency-dependent Ra- 
man cross section, see Fig. \S\ (b), and the frequency- 
dependence of the electromagnetic enhancement as 
shown in Fig. |21 The free-molecule Raman cross sec- 
tion is considerably higher for both hflL = 2.4 eV and 
TI^Il — 2.5 eV than for hrtL — 2.6 eV, which explains 
why the Raman peak in the latter case is so small. 
In addition, as a result of the resonant maximum at 
« 2.2 eV, in Fig. |2 the combined EM enhancement 
|M(Ol - Wvib)P|Af(f]L)P is the largest for = 2.4 
eV and smallest for hil^ = 2.6 eV. 

Figure El shows the spectra that result when the de- 
phasing rate is varied. The diagram only shows the por- 
tion of the spectrum where the Raman peaks (the prin- 
cipal one and the first overtone) falls. Both these cross 
sections are, at least with a laser photon energy as close 
to resonance as 2.45 eV, quite sensitive to the dephas- 
ing rate. This is particularly apparent for the Raman 
overtone peak. For creation of multiple vibrational exci- 
tations to take place, the molecule must have a chance 



FIG. 9: (color online) Calculated spectra for four different 
dephasing rates. The molecule-particle separation is d = 5 
A, and the rest of the parameter values are the same as in 
Fig. |7| The curves have been shifted vertically for increased 
clarity. The truncated peak for 'jph = 50 x 10^^ s"'^ has height 
^ 200 X 10-i''cmV(meVsr). 
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FIG. 10: Calculated spectra when the molecule is placed 
near only one Ag particle with radius 400 A. The molecule 
parameter values are the same as in Fig. [7| 



to "investigate" the displaced oscillator potential during 
a longer time than in the case of single-excitation cre- 
ation. This means that the multiphonon processes are 
more easily disrupted when the dephasing rate goes up. 

In Fig. ^1 we show results for a different nanoparti- 
cle configuration than the one considered so far. The 
results here refer to a molecule placed near only one 
silver particle. In this case there is some electromag- 
netic enhancement M, because the metal particle can 
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be (resonantly) polarized. But due to the lack of elec- 
tromagnetic particle-particle coupling the enhancement 
has a much smaller magnitude than in the two-particle 
case. Typically |Mp reaches values of 20-40 (depend- 
ing on frequency) for c? = 5 A. M is less sensitive to d 
than in the two-sphere case, the most important contri- 
bution coming from the dipole field around the particle, 
i.e. \M\ ~ /{R + dY . However, the damping rates 
due to losses in the particles, especially electron-hole pair 
creation, are still comparable to those in the two-particle 
case. This means that fluorescence is very strongly sup- 
pressed for small d when there is only one Ag particle. 
For example, for d = 5 A a rough estimate of the ratio 
|M|4/|Mdp with |M|2 « 30 and jM^P taken from Fig.EI 
shows that we can expect a suppression (or quenching) of 
the fluorescence by a factor of ^ 10'^ compared with the 
free- molecule case. In Fig. ^]we see in fact a difference 
by about 3 orders of magnitude between the fluorescence 
when c? is 5 A and when d is 1000 A, the latter d yields 
situation quite similar to the free-space case. For d — ^ 
A the fluorescence is suppressed to the extent that Ra- 
man peaks stand out from the background, however, the 
absolute Raman cross sections are of course too small to 
be detectable in an experiment with a single molecule. 

With increasing d in Fig. ^| the fluorescence yield 
first increases because \Md\ decreases rapidly whereas 
\M\ falls off at a much slower rate. Then when d reaches 
values of 50-100 A ( see the 70-A curve) the fluorescence 
has a maximum, because the rate of decrease in \M\'^ 
overtakes that of iM^p. For even larger d jMp « jM^p, 
so the fluorescence cross section behaves as |Mp there. 
The spectrum calculated for d = 1000 A approaches what 
one finds for the molecule in free space. 

The results presented so far have been calculated with 
moderate laser intensities. We have used the value 
10'*V/m for in Eq. which corresponds to /in ~ 
0.13 ^W/(/im^). This intensity is so low that, even in 
spite of the EM enhancement, the molecule spends al- 
most all the time in the ground state (assuming, as we 
have done before, that thermal excitations can be ne- 
glected). The results presented in Fig. 1111 fwhere we re- 
turn to a situation with two nanoparticles) have been 
calculated with considerably higher incident intensities. 
This means that we encounter situations where there is 
an appreciable probability of finding the molecule elec- 
tronically and/or vibrationally excited. For very strong 
incident fields also the response properties of the molecule 
changes, essentially because the probability of finding 
the system in a certain level becomes time-dependent 
(as a result of Rabi oscillations) . In our calculations the 
stationary-state density matrices pss a-nd po describe a 
time-average, the way the molecule "looks" on the aver- 
age after long time. The dipole-dipole correlation func- 
tion (p'-^-' (0)p'-+^ (i)), on the other hand, contains infor- 
mation about the molecule dynamics over a shorter pe- 
riod of time, starting from a certain initial state. Here 
Rabi oscillations due to excitation and deexcitation by 
a strong driving, external field show up in the results. 




Photon energy (10 meV/div.) 



FIG. 11: (color online) Calculated spectra for a molecule 
placed between two nanoparticles, with d = 5 A, for a number 
of different incident intensities. The other parameter values 
are the same as in Fig. 13 Only the energy ranges around 
the Raman Stokes peak at Tiu)s = 2.29 eV and the Raman 
anti-Stokes peak at Tiujas = 2.61 eV are shown. 



along with other, more apparent, time-dependent aspects 
of the molecular dynamics such as dipole oscillations and 
damping. 

In Fig. we restrict the attention to parts of the 
spectra in the frequency ranges around the Stokes and 
anti-Stokes peaks with center frequencies ujs = ^ge—^v\h 
and uJAs = f^ge+Wvib, respectively. The different spectra 
correspond to intensities ranging from w 13/iW/(/xm^) 
(at Eq = lO^V/m) to « 3.0mW/(/im2) (at Eq = 
15 X lO^V/m). The rest of the parameter values have 
been chosen the same way as in the calculation behind 
Fig. [7| At the lowest intensity we get a Stokes peak that 
is nearly identical to the one in Fig.EI but with increasing 
intensity the peak height diminishes quite rapidly, and it 
is also broadened. At the same time, at the anti-Stokes 
frequency, a small bump for Eq = 10^ V/m develops into 
a marked peak at £^o = 5 x 10^ V/m, and subsequently 
also this peak is reduced in height and broadened with 
increasing Eq. Thus, with incident intensities of the or- 
der of 1 mW//im^ it is possible to obtain anti-Stokes Ra- 
man scattering (in the model) with a cross section that 
is experimentally detectable. The peak values found in 
the figure must be multiplied by an effective peak width 
~ lOmeV and the effective solid angle (for dipole scat- 
tering) 87r/3 to give the total Raman anti Stokes cross 
section. The fact that the anti-Stokes signal is weaker 
than the Stokes signal is here partly due to the choice of 
laser frequency which means that the Stokes peak falls 
at a maximum in the EM enhancement while the anti- 
Stokes peak is at a minimum for |Mp. 

To discuss the tendencies in Fig. ^] as a function of 
the incident field it is useful to calculate the probabili- 
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FIG. 12: The probability of finding the molecule in differ- 
ent states and for a number of different driving fields. States 
numbered through 5 on the horizontal axis are in the elec- 
tronic ground state, i.e. |(/)0; n) for < n < 5, whereas states 
6 through 11 correspond to the electronically excited states 
|e)a::o;0) through |e)a;o;5). 



ties of finding the molecule in different states. This is 
shown in Fig. 1121 The development at moderate inten- 
sities can be well understood just from looking at the 
probabilities of finding the molecule in either of the two 
lowest states |(7)|0;0) and |(ji)|0;l). With increasing in- 
tensity it becomes more likely to find the molecule in the 
vibrationally excited state which is a good initial state 
for the anti-Stokes Raman process. At the same time 
the vibrational ground state, the most common initial 
state for a Stokes Raman process, becomes more and 
more depopulated leading to a decreased intensity there. 
However, comparing the probabilities in Fig. 1121 with the 
changes in peak heights in Fig. ^2 it becomes clear that 
the behavior at higher intensities cannot be understood 
exclusively in terms of the probabilities. If this was true, 
the Stokes peak would not diminish as rapidly as it does, 
and the anti-Stokes peak would continue to increase in 
height. The reduced strength of both the anti-Stokes and 
Stokes peaks is instead the result of the molecule response 
becoming time-dependent at strong driving fields which 
leads to a smearing-out of all spectral features. 



tions of order 10~^^cm^, thus comparable to the values 
found in single molecule SERS experiments provided 
the laser frequency is reasonably close to an electronic 
excitation frequency of the molecule, and the electromag- 
netic environment (the nanoparticles) provide an EM en- 
hancement of the Raman cross section by some 10 orders 
of magnitude. Thus, we conclude that the 14 to 15 orders 
of magnitude enhancement sometimes cited in connec- 
tion with single-molecule SERS may neither be present 
nor needed for the effect to occur. 

Our study gives an opportunity to study how the Ra- 
man and fluorescence parts of the spectra develop as 
the electromagnetic enhancement is varied. Most of the 
trends found here agree well with what one can expect 
compared with experiment. For the highest enhance- 
ments the Raman peaks stand out from the fluorescence 
background. The rather rapid decrease of the Raman 
cross section with increasing molecule-particle separation 
agrees qualitatively with results found in experiments 
where the metal particles are coated with a thin layer 
of dielectric material, the thickness of which can be con- 
trolled with reasonable precision At the same time 
the fluorescence background is more structured here than 
in most experimental results found in the literature^*^ 
although some exceptions do exist 4^ 

The quenching of the fluorescence (Fig. I10|) once a 
molecule placed near a metal surface is a well-known phe- 
nomenon, and systematic studies have also found a local 
maximum in the fluorescence cross section as a function 
of the molecule-surface distancCfM Compared with this 
study, the fluorescence maximum occurs closer to the sur- 
face in our results because of the faster drop of the elec- 
tromagnetic enhancement with increasing distance from 
a finite nanoparticle than from a flat surface. 

The model we have used in this work is a basic one con- 
taining a minimum of ingredients. To better account for 
some aspects of the SERS phenomenon, further develop- 
ments are needed: (i) More vibrational modes should be 
added to better describe a fairly large organic molecule, 
(ii) More electronic levels could also be included. This 
could be a way to model "blinking" phenomena due to 
the molecule's spending some time in a metastable state 
that is not directly optically active. One could also in- 
troduce electronic states leading to dissociation of the 
molecule, thus modeling photobleaching effects, (iii) Ul- 
timately one should also try to include electron transfer 
processes between the molecule and the metal particles. 



VII. DISCUSSION AND SUMMARY 
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We have performed a theoretical analysis of surface- 
enhanced Raman scattering and fluorescence based on 
a treatment that combines electromagnetic enhance- 
ment effects and molecule dynamics. We used a simple 
molecule model, however, with parameter values chosen 
in a realistic way. 

The calculated results give Raman scattering cross sec- 
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APPENDIX A: ELECTROMAGNETIC 
CALCULATION 

The coefHcients c in Eq. H21|l are explicitly given by 

1 

(Al) 

when T = t' , and 

g 2i' + 1 

■ ' ■ ' ' i^e'if + 1) ^.f-^ 

'^^hn ^i.rn+j4'-l.m'+j{R) C'l'',„i',f'-l,„i'+j i (A2) 

when T = t' . In these expressions 

^ and - Vi^Tm)ii±m + l) 



F 



(A3) 



7^+l = -v/47r/3, while 7^o = TZ-i = V47r/3, (A4) 
and the Gaunt integrals Cf^^, ,„/ are defined by 

Q":™/',™' = / dnYi^{n)YLMmyi'm'{^). (as) 

The coefficients A in Eq. IjAip are given by 

i+f.' 

Ae,mr.rn'{R) = 4vr ^ i-iy-"^-'' hL{k\R\) 

L = \l-l'\ 

^Yl,^^^,{R)CI:Z7,%,. (A6) 

Finally, the coefficients describing the external field in 
Eq. are 



'lfm,Gxt 



47r/e*^~= ^ r,%^./fc)F/„[e*.^o], (A7) 



and 



(A8) 



where the unit vectors e_i, eo, and ei are given by 

eo = 5 and e±i = {x±y/i)/V2, (A9) 

and /c = A:/|fc| is a unit vector in the direction of k. When 
instead the spheres are driven by a dipole zpoe~'"' at the 
molecule position we have 



'^T^m, oxt 



(AlO) 



APPENDIX B: ELECTRON-HOLE PAIR 
DAMPING 

In this Appendix we outline the calculation of the 
electron- hole pair contribution Peh/Pfree to the damp- 
ing rate enhancement |Mrfp. The non-local dielectric 
response of the metal particles is treated within d- 
parameter theory^^^ and the derivation to a large extent 
follows Ref. 0. 

The idea is to model the interaction between the molec- 
ular dipole and the degrees of freedom (electron-hole 
pairs) of the nearby metal particles by a linear coupling 
between the dipole and a number of boson modes. To be 
specific, the dipole points along the z axis and is placed 
at z = h between two fiat metal surfaces at z = and 
z = L. The interaction Hamiltonian can then be written 



Flint ^ ^\ Cq^(yhq^QL 



(Bl) 



Here p is the molecule dipole operator, 6 and are anni- 
hilation and creation operators for the bosons (with in- 
plane wave vector q and another branch index a) and the 
coupling coefficients C^,a are dependent on the position 
h of the dipole. Using the Fermi golden rule, Hint now 
gives a decay rate from the excited state to the ground 
state of the molecule which can be written 

w^^\p,\''A I ^Y.\c^,^\H{n^i,,-nu^^^^). (B2) 

The above expression is only useful if we have a way of 
calculating the coefficients C^.q. To do this we evaluate 
the energy dissipation from a classical dipole placed at 
the position of the molecule to the bosonic degrees of free- 
dom. Again, this expression will contain the coefficients 
Cq-^ct, but with a classical dipole the energy dissipation 
rate can also be calculated using classical electrodynam- 
ics, and expressed in terms of geometric parameters and 
the dielectric function of the metal. We write the classi- 
cal dipole moment as 



(B3) 



By letting p(t) take the place of p in the Hamiltonian, 
we get, from the Fermi golden rule, an energy dissipation 
rate 

W = 2T:uj\pi\'^^\Cij,a?5{huj - huJij^a)- (B4) 

Next we must calculate the dissipated power W within 
the framework of classical electrodynamics. To this end 
we place the dipole with a dipole moment given by Eq. 
(|B3|) between the metal surfaces, dX z = h. Since all the 
distances involved in this calculation are very short we 
can ignore effects of retardation and express the solution 
in terms of a scalar potential (j){t) = 0ie~*'^* -I- (^*e''^*, 
where 4>i can be expressed as a Fourier transform 



(j>i{x,y,z) 



(2^) 



0,(z)e^'^-'-|i, 



(B5) 
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(r|| = XX + yy) . For <f> to satisfy Poisson's equation 

V'^cf = -p/eo, (B6) 

(where p is the charge density due to the dipole), (t>q{z) 
must be a hnear combination of two exponentials e^n ^ and 
e~'ii^. In the classical calculation we treat the surfaces in 
terms of their surface response functions: gi(q,u}) for the 
surface at 2 = and g2{q,iLj) for the surface aX z — L. 
The surface response function gives the negative ratio 
between the "reflected" potential (decaying as one leaves 
the surface) and the "incident" potential (decaying as 
one approaches the surface. Thus, 



(B7) 



where the potentials should be evaluated right at the 
surfaces. Within d-parameter theory the surface response 
function is given to first order in qii = a: 



^.32 



9(9,'^) = 



e{Lj) + 1 



e(c.) 



(B8) 



The local dielectric function £{lo) is the same as used in 
the Mie calculation and d^{uj) is the d-parameter func- 
tion. In our calculations we have evaluated Im[(ij_] from 
Table I of Ref. 32 using rg = 3 and Tiujp = 9 eV appro- 
priate for silver in the low-frequency regime. The real 
part Re[(i_L] plays a less important role; we set it to the 
constant value Re[(ij_] = 1 A here. 

The solution for (pi can now be expressed in terms of 



either the total incident potential at z = 0, 0g.down, as 



<j)i{x,y,z) 



'<t>qAo^n (e«" - 5ie-«") , (B9) 



or, by a similar expression, in terms of ^g.up, the incident 
potential at the upper interface. 0q,down and 0q,up are 
given by 



0(j,down 



and 



Pi 



-q(L + h) 



2eo 1 - ffiff2e" 



-2qL 



(BIO) 



(Bll) 



The Poynting vector S = E x H at the two interfaces 
can here be approximated by an expression involving only 
(f>; the dissipated power is then given by 



W = 4e, 



(fq 



rujq 



|09,downplmgi 



.,upplmg2 



(B12) 



But this power should equal the one found in Eq. (|B4p . 
and in this way we get a relation between the classical 
quantities found here and the sum over a in Eq. (|B4|I . 
(the sum over q goes over to an integral). Inserting the 
so obtained expression into Eq. I|B2|I gives 



,2 r 



dq\\ 2 
(2^^°'^ll 



2'/>q,dc 



Pi 



Imgi 



r 



pi 



Img2 



(B13) 



The electron- hole contribution to \Md\'^ in Eq. H26|) can 
now be found as 



Peh/Pircc = w/7rad,0, 



(B14) 



where 7rad,o is found from Eq. Ml\ with |/p = 1. 

The integrand of Eq. I|B13|I has been plotted, for two 
different photon energies, as a function of q in Fig. 1131 
both with and without a nonzero imaginary part for the 
function di_{uj). This illustrates different contributions 
to the damping rate. At small q we have damping medi- 
ated by relatively long-wavelength interactions between 
the molecule and the substrates. This part of the damp- 
ing is actually captured by the Mie calculation discussed 
in the main text. In this range of q space we see a fairly 
sharp peak around 0.5 nm~^ that is due to losses to 
a resonant interface plasmon mode formed between the 
two metal surfaces. We also see that the use of a non- 
local dielectric function, a nonzero value for Imdj^(ci;), 



has essentially no effect here. Moving towards higher q 
values we encounter contributions to the damping that 
are not included in the Mie calculation. There are two 
reasons for this: (i) It becomes technically difficult to 
go to very high values for the angular momentum i and 
consequently rapid (high q) variations of the fields are 
not accounted for. (ii) At these larger wave vectors the 
non-local effects (Imd^(w) 7^ 0) that are not so easily 
included in the Mie calculation provides the dominant 
contribution to the damping. This is clearly seen in Fig. 

In order not to double-count contributions to the 
damping already included in the Mie calculation we em- 
ploy a low-g cutoff in Eq. (|B13|) using gmin = Inm^^. 
It is not really possibly to determine an exact position 
for the cutoff, since there is no exact correspondence be- 
tween angular momenta and wave vectors, but 1 nm~^ 
is a reasonable value to use together with i? = 40 nm and 
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